Oxygen impurities in NiAl: relaxation effects 
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We have used a full-potential linear muffin-tin orbital method to calculate the effects of oxygen 
impurities on the electronic structure of NiAl. Using the supercell method with a f 6-atom supercell 
we have investigated the cases where an oxygen atom is substitutional^ placed at either a nickel 
or an aluminum site. Full relaxation of the atoms within the supercell was allowed. We found that 
oxygen prefers to occupy a nickel site over an aluminum site with a site selection energy of 138 mRy 
(21,370 K). An oxygen atom placed at an aluminum site is found to cause a substantial relaxation 
of its nickel neighbors away from it. In contrast, this steric repulsion is hardly present when the 
oxygen atom occupies the nickel site and is surrounded by aluminum neighbors. We comment on 
the possible relation of this effect to the pesting degradation phenomenon (essentially spontaneous 
disintegration in air) in nickel aluminides. 
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Nickel aluminides are potentially important industrial 
alloys in high-temperature structural applications, such 
as gas turbine engines, because they possess many desired 
physical properties. Among these are their low density, 
high strength, and good thermal conductivity, the last 
being an important factor in dissipating the heat gener- 
ated by the engines. The great interest in these alloys 
has resulted in a large body of literature. jlr|| 

Many properties of NiAl have been intensely investi- 
gated using various computational electronic structure 
methods. The band structure of NiAl has been experi- 
mentally studied using photocmission and calculated us- 
ing the linear augmented Slater-type orbital method by 
Lui et al. |7j They found that NiAl behaves like a good 
itinerant metal: the self-energy corrections in NiAl are 
significantly less than in pure ferromagnetic nickel. Band 
structure results should therefore provide a good descrip- 
tion of the electronic properties of NiAl. Kim et al. 
H also calculated the band structure of NiAl using a 
semirelativistic linearized augmented plane wave method 
and compared the result with the experimental optical 
spectra. [|]|| Their result agreed with the result of Lui, 
and they found no need to incorporate self-energy cor- 
rections into the spectrum in order to fit their experi- 
mental data. This is markedly different from the result 
for CoAl where the corrections are needed to make a 
good fit. M One can therefore put some confidence in 
the band-structure results for NiAl. 

Most research on nickel aluminides has been devoted to 
the effort to understand the complex effects of impurities 
and other kinds of crystal defect, pi Uses of nickel alu- 
minides, especially in polycrystalline form, as structural 
materials are limited by their brittleness at high temper- 
ature (over 800 °C) where they become prone to brittle 
intergranular fracture and by their low ductility at inter- 
mediate temperatures. Certain additive impurities have 
been found to greatly improve the room-temperature in- 
tergranular cohesion and the tensile ductility. The most 
notable of these cohesion enhancers is boron fill] which 
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FIG. 1. The 16-atom supercell used in the present work. 
The cubic supercell is built up from 2 3 unit cells of NiAl. The 
impurity atom (O) is placed at the center of the supercell, 
substituting for Ni in this picture (supercell NirAlsO). We 
also use a similar supercell, where the oxygen atom is placed 
at an Al site (supercell NisAlrO). 



can improve the tensile ductility by an order of magni- 
tude and change the fracture mode from intergranular 
to transgranular. 10 Sun et al. |12) have studied the ef- 
fects of boron and hydrogen on the cohesion of M3AI us- 
ing a full-potential linear muffin-tin orbital (FP-LMTO) 
method. They concluded that boron improves the local 
cohesion by reducing the bonding charge directionality 
around the nickel atoms and by inducing an increase of 
the interstitial bonding charge. 

At intermediate to high temperature, nickel alu- 
minides, and many other intermetallics, exhibit brittle 
intergranular fracture due to the oxygen-induced embrit- 
tlement. |l^| At high temperature, the oxidation will se- 
lectively attack the least noble constituent, which is alu- 
minum, and form the oxide product AI2O3. CjJ In some 
cases, these oxide products can provide a stable oxide 
layer that protects the alloy underneath it from further 
oxygen attack. NiAl is one of the intermetallics that are 
able to form a protective oxide layer and therefore is used 
as a coating for other intermetallic alloys. p3[ 



TABLE I. Equilibrium lattice constant a r for fully relaxed 
supercell, its total energy E r , and the corresponding relax- 
ation energy Er (relative to the unrelaxed supercell). 



Supercell 



(a.u.) 



E r (Ry) 



E R (mRy) 



Ni 8 Al 7 
Ni 7 Al 8 



10.5878 
10.6567 



-27,831.7387 
-25,278.8899 



-0.675 
-1.236 



The oxidation resistance decreases as we move toward 
the nickel-rich part of the Ni-Al phase diagram. The 
most damaging effect of oxygen occurs when it causes 
the pesting degradation phenomenon (essentially spon- 
taneous disintegration in air). Some intermetallic com- 
pounds that form protective coatings at high tempera- 
ture literally disintegrate when heated in the intermedi- 
ate temperature range. [ |13|Jl4| This occurs only on poly- 
crystalline samples, so it is sometimes also termed in- 
tergranular attack. C3] Among nickel aluminides, pest- 
ing phenomenon has been observed in NiAl pUla] and 
Ni 3 Al. jT§ 

Despite the extreme effects that it can cause, we are 
unaware of any previous ab-initio study performed on the 
effects of oxygen in nickel aluminides. Moreover, most 
studies of other impurities are limited by not allowing the 
atoms within the supercell to relax. Relaxation effects 
are especially relevant to the case of impurity atoms lo- 
cated at the grain boundary where they can be relatively 
free to move and bond to certain constituent atoms of 
the host alloy. 

In this work we have used a full-potential linear muffin- 
tin orbital method within the local density approxima- 
tion (LDA) to study the effects of oxygen impurities 
on the electronic structure of NiAl. The details of this 
method have been documented elsewhere. fll7|-|l9|] For the 
work reported in this paper, we use nine (spd) orbitals 
for each atom in the supercell and assign three LMTO 
basis functions for each orbital, corresponding to k — 
— 1.50,-0.50, and 1.00, respectively. Here k is the wave- 
number parameter; the square of k is the absolute value 
of the energy of the basis function, measured relative to 
the muffin-tin zero, while its sign is equal to the sign of 
the energy. [EOJ The large number of basis functions used 
provides our full-potential method with good flexibility 
to find the lowest-energy density. 

The 16-atom supercell that was used in the computa- 
tion is shown in Fig. EJ. Pure NiAl crystallizes in the B2 
structure with lattice constant of 5.4450 a.u. {2lJ The cal- 
culated equilibrium lattice constant using our FP-LMTO 
method is 5.3451 a.u which agrees with the experimental 
value within 2%. The supercell is constructed from 2 3 
unit cells of NiAl. The oxygen impurity atom is placed 
at the center of the supercell, replacing a nickel or an 
aluminum atom. For a given supercell lattice constant, 
we allow the atoms in the unit cell to relax to find the 
minimum total energy for that lattice constant. By sym- 



metry, only the eight nearest neighbors of the oxygen 
atom are allowed to relax, and they can only move radi- 
ally away or toward the center oxygen atom (we do not 
consider the possibility of symmetry breaking). Only one 
parameter is needed to describe the relaxation, namely 
the distance of these neighboring atoms from the central 
oxygen atom. 

The computed total energies for the relaxed 16-atom 
supercell are shown in Fig. and compared with the 
curves for the unrelaxed supercell. Each continuous 
curve shown in Fig. pi is obtained by fitting 8-9 LMTO 
data points computed around the minimum-energy lat- 
tice constant using an eight-parameter fitting function. 
The calculated relaxation energies Er are listed in Ta- 
ble H. We use the following definition for the relaxation 
energy 



Er — E r (a r ) — E u (a u ), 



(1) 



with a u and a r being the lattice constant that give the 
minimum energy for the unrelaxed (E u ) and relaxed 
(E r ) configurations, respectively. The computed energy 
curves lie very close to each other, and the relaxation 
energies are relatively small, being about 100 K and 200 
K (in equivalent temperature scale) for oxygen at the Al 
and Ni sites, respectively. These numbers are about the 
same as the accuracy of the LMTO-LDA method, which 
has been estimated at about 1 mRy (158 K). j22| Numer- 
ically, however, our FP-LMTO method did not produce 
any notable fluctuation and we believe that some confi- 
dence can be placed in these numbers. 

The computed relaxation data are shown in Fig. 
where we track the nearest-neighbor distances between 
the atoms in the unit cell on the (011) plane as we change 
the supercell lattice constant. Figure 0(b) displays the 
results for the case of oxygen at the Ni site. Shown are 
d(Al— O), which is the distance between the central O 
atom and one of its eight nearest-neighbor Al atoms, and 
d(Ni— Al), which is the distance between the Al atom and 
its nearest neighbor Ni atom (at one of the corners of the 
supercell). Without relaxation, the oxygen atom will be 
at (0,0,0), the aluminum atom at, e.g., (j,j,j), and 
the nickel atom at (~, |, |). As we vary the lattice con- 
stant a, the distances scale linearly with it: d(Al— O) = 
d(Ni— Al) = a\/3/4. This relation is still approximately 
followed when we relax the atoms, as shown in Fig. 0, 
with the aluminum atoms only attracted slightly more 
toward the central oxygen atom. 

A significantly different situation occurs if we place 
the oxygen atom at the Al site, as shown in Fig. 0(a). 
In this case, the nearest neighbors of the central O atom 
are Ni atoms, the distance between them is d(Ni— O), 
and each of these Ni atoms neighbors an Al atom at its 
nearest corner of the supercell, where the distance be- 
tween them is denoted by d(Ni— Al). Under pressure, for 
lattice constants smaller than the equilibrium value, the 
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FIG. 2. The FP-LMTO total energy with 16-atom super- 
cell for the case of oxygen substituting for an aluminum atom 
(a), and oxygen substituting for a nickel atom (b). In both 
panels, the top curve is for an unrelaxed supercell while the 
bottom curve is for the relaxed supercell. The equilibrium 
lattice spacing and the relaxation energy are listed in Table 
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FIG. 3. Nearest-neighbor (n.n.) distances between atoms 
in the unit cell on the (011) plane as a function of the supercell 
lattice constant, (a) Oxygen at Al site. For small lattice 
constants the movement of the Ni atoms is toward the central 
oxygen atom and away from the corner Al atoms. At large 
lattice constants the direction of movement reverses: the Ni 
atoms are displaced significantly closer to the corner Al atoms. 
This presumably reflects a steric repulsion between oxygen 
and nickel in the NiAl environment, (b) Oxygen at Ni site. 
Apparently, the Al atoms are attracted about equally strongly 
by the O and the Ni atoms. This results in small relaxation. 



Ni atoms slightly relax toward, the central oxygen atom, 
away from the corner Al atoms. However, as the lattice 
constant is increased above the equilibrium value, the in- 
tervening Ni atoms very quickly start to move away from 
the central oxygen atom and relax closer to their neigh- 
boring Al atoms. We view this as reflecting a mutual 
repulsion between the oxygen and nickel atoms in the 
NiAl environment. Experimentally, the formation rate 
of aluminum oxide during exposure of nickel aluminides 
to oxygen is known to be much higher than that of nickel 
oxide. ^|JI|] 

The relaxation behavior of the atoms is thus seen to 
differ markedly depending on whether the impurity oxy- 
gen atom occupies an Al or a Ni site, especially in the 
stretched-supercell case, where the atoms have greater 
freedom to move around. This may have some relevance 
to the case of oxygen attack on polycrystalline nickel 
aluminides (pesting). We believe this reflects a situa- 
tion where the oxygen atoms cause a sort of "wedge ef- 
fect", seeping into the grain boundaries to reduce the 
intergranular cohesion and opening up the polycrystal 
wider for more infiltration of oxygen. This provides a 
self-propelling mechanism for oxygen atoms to infiltrate 
a polycrystal of nickel aluminide and, along the way, de- 
stroy the intergranular cohesion between the grains, effec- 
tively disintegrating the polycrystal. In this scenario the 
pesting phenomenon can be seen to be fueled by the com- 



bination of two major factors: the thermal intergranular 
diffusion of oxygen and the strongly preferential bond- 
ing of oxygen with one of the components of the alloy 
(aluminum in nickel aluminides). More insights into this 
phenomenon can perhaps be obtained by performing a 
molecular-dynamics (MD) simulation using microscopic 
parameters that are extracted from an ab-initio calcula- 
tion such as reported in this paper. Campbell et al. have 
performed such an MD simulation for the oxidation of 
aluminum nanoclusters. p3| 

The energetics of the preferential bonding of oxygen 
in NiAl can be studied by calculating the site selection 
energy. [^4J The total energy for pure bulk metal (per 
atom) is used to provide a reference energy for the con- 
stituent species. In our case, these are the total energies 
of fee Ni and fee Al. These have been calculated using 
the same FP-LMTO method and their values are 



£(Ni) = -3, 036.8304 Ry, 



E(M) = -483.8440 Ry. 



(2) 



(3) 



The total energies for the 16-atom supercell with one 
oxygen atom replacing a Ni or an Al atom are listed in 
Table EL We list the minimum of the total energy of the 
relaxed supercell. The site selection energy is defined to 
be the difference between the following two values p3] : 
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Fermi Energy: 
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FIG. 4. Total density of states for the 16-atom supercell 
calculated using the FP-LMTO method. For comparison, in 
each panel we also show the negative of the density of states 
for pure NiAl. Vertical arrows on the energy axis point to 
the positions of the Fermi energy, (a) Oxygen at Al site, (b) 
Oxygen at Ni site. 



E m = E(Ni) + £(Ni 7 Al 8 0) = -28, 315.7204 Ry, (4) 




Fermi Energy: 

NiAl : 0.8768 
NLALO: 0.8613 



DOS is the existence of the sharp peaks due to the d- 
orbitals of nickel, which hybridize only weakly with the 
other orbitals. Placing the oxygen atom at the Al site 
in NigAlyO allows for some hybridization between these 
d-orbitals and the delocalized p-orbitals of oxygen. This 
results in reduced sharpness of the Ni d-state peaks in 
the DOS without essentially any shift in the position of 
the peaks, as can be seen in Fig. 0(a). On the other 
hand, if we place the oxygen impurity at the Ni site, 
then the nickel atoms will not have the oxygen atom as 
their nearest neighbor. The oxygen atom will have the 
aluminum atoms as its nearest neighbors. Being the more 
electronegative element, oxygen will interact with the va- 
lence electrons of Al and will localize a portion of those 
electrons around itself. This depletes the Ni sites of some 
of the electrons from the Al that was formerly occupy- 
ing its site in the pure NiAl case. The end result of this 
is a lower electrostatic potential at the Ni sites and the 
reduction of the onsite energies of the (i-orbitals of Ni, 
without much change in their spatial extent. This trans- 
lates to an almost rigid downward shift in the position of 
the Ni d-state peaks in the DOS without much alteration 
in their width. This downward shift of about 25 mRy can 
quite readily be discerned in Fig. 0(b). 
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E M = E(M) + £(Ni 8 Al 7 0) = -28, 315.5827 Ry. (5) 

Taking the difference, an oxygen atom in NiAl will prefer 
to occupy a nickel site with the site selection energy of 



AE = 137.7 mRy = 21, 730 K. 



(6) 



That oxygen will prefer to occupy a nickel site over 
an aluminum site seems to contradict one's expectation 
based on the atomic radii of the constituent atoms. Since 
the atomic radius of oxygen is closer to that of aluminum 
rather than nickel, one would expect that the oxygen 
would prefer to occupy an aluminum site over the nickel 
site. To understand our result, we need to recall the 
previously deduced mutual repulsion between oxygen and 
nickel in the NiAl environment. An oxygen atom will 
therefore prefer to be surrounded by nearest neighbors of 
aluminum, rather than nickel, atoms. It will achieve this 
simply by occupying a nickel site. 

Finally, in Fig. we show the total density of states 
(DOS) for the two supercell systems that we study in 
this work, NigA^O and NirAlgO, calculated from the 
FP-LMTO energy bands using the tetrahedron method. 
For comparison, we also show the corresponding total 
DOS for pure NiAl. The dominant feature of the NiAl 
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